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ABSTRACT 

We present results from the first robust investigation of galaxy morphology as a function of envi- 
ronment at z > 1.5. Our study is motivated by the fact that star- forming galaxies contained within a 
protocluster dX z — 2.3 in the IIS1700+64 field are characterized by significantly older ages and larger 
stellar masses on average than those at similar redshifts but more typical environmental densities. In 
the analysis of HST/ ACS images, we apply non-parametric statistics to characterize the rest- frame 
UV morphologies of a sample of 85 UV-selected star- forming galaxies at z = 1.7 — 2.9, 22 of which 
are contained in the protocluster. The remaining 63 galaxies are not in the protocluster but have 
a similar mean redshift of z 2.3, and constitute our control sample. We find no environmental 
dependence for the distributions of morphological properties. Combining the measured morphologies 
with the results of population synthesis modeling, we find only weak correlations, if any, between 
morphological properties and measures of stellar population such as stellar mass, age, extinction and 
star-formation rate. These findings are similar to results from other recent work. Given the incomplete 
census of the protocluster galaxy population, which is weighted towards star-forming galaxies, and 
the lack of correlation between rest-frame UV morphology and star-formation history at z ~ 2 within 
our sample, the absence of environmental trends in the distribution of morphological properties is not 
surprising. We do, however, find some non-environmental trends with morphology. The degree of 
nebulosity and apparent rest-frame UV size appears to correlate with the UV-luminosity (uncorrected 
for dust extinction), which we demonstrate is not an artifact of finite signal-to- noise. Additionally, 
using a larger sample of photometric candidates, most of which do not have spectroscopic redshifts, 
we compare morphological distributions for 282 UV-selected and 43 near-IR-selected galaxies. While 
the observed difference in the degree of nebulosity between the two samples appears to be a byproduct 
of the fainter average rest-frame UV surface brightness of the near-IR-selected galaxies, we find that, 
among the lowest surface brightness galaxies, the near-IR-selected objects have significantly smaller 
angular sizes. 



Subject headings: galaxies: high redshift 
observations 



galaxies: structure — galaxies: clusters — cosmology: 



1. INTRODUCTION 

Galaxies populate the local universe in a very regular 
way. Dynamically hot elliptical galaxies are found prefer- 
entially in large galaxy clusters, especially the very cen- 
ters (Dressier 1980; Postman & Geller 1984; Goto et al. 
2003). Dynamically cold spiral galaxies tend to lie in 
the less dense environments of groups or the field. The 
name given for these trends is the "morphology-density" 



relation. The morphology-density relation is already in 
place by z = 1, although there are quantitative, but not 
qualitative, differences between z = 1 and z — (Post- 
man et al. 2005; Smith et al. 2005). For example, the 
elliptical fraction in clusters is lower at z = 1 than at 
z = 0, but the cluster elliptical fraction is higher than 
the elliptical fraction in the field in both epochs. One of 
the outstanding problems in the theory of galaxy evolu- 
tion and structure formation is understanding when and 
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how this relation is established. 

The morphology-density relation is entangled with sev- 
eral other observed relations between physical properties 
of galaxies and their classification and environment. In 
the local universe, elliptical galaxies are structurally dif- 
ferent and, on average, have higher mass, less active star 
formation (as well as less gas available for star forma- 
tion), and redder colors than spiral galaxies (Kauffmann 
et al. 2003). Therefore, in addition to a morphology- 
density relation, strong trends of galaxy color with en- 
vironment are observed in the Sloan Digital Sky Sur- 
vey (SDSS; Kauffmann et al. 2004; Blanton et al. 2005). 
These trends have been traced out to z = 1.3 with the 
DEEP2 galaxy redshift survey (Davis et al. 2003; Cooper 
et al. 2006b,a; Gerke et al. 2006). Studies of the local 
universe strongly imply that the correlations between 
color and star formation history and environment are 
more fundamental than the morphology-density relation 
(Kauffmann et al. 2004; Blanton et al. 2005). As part 
of the program of understanding galaxy evolution, it is 
important to explain the relation between the physical 
structure of galaxies and their star formation history, and 
the correlations between the various galactic properties 
and environment. 

There are many theories to explain how such correla- 
tions arise as part of the hierarchical picture of struc- 
ture formation. Several authors have considered how 
much of the morphology-density relation is simply a re- 
sult of the merger histories of the host halos. Benson 
et al. (2001) graft semi-analytic models onto a cosmo- 
logical dark matter N-body simulation (Jenkins et al. 
1998; Kauffmann et al. 1999) and examine the result- 
ing morphology-density and color-density relations at the 
current epoch. They find qualitative agreement with ob- 
servations of the early-type fraction and average galaxy 
color as a function of environment. Their semi-analytic 
model overpredicts the red fraction at the center of clus- 
ters, and underpredicts the number of spheroids in the 
field at 2 = 1 (Benson et al. 2002). Okamoto & Na- 
gashima (2001) also use semi-analytic models to inter- 
pret their dark matter simulations, but focus only on 
the morphology-density relation within clusters. While 
they can reproduce the elliptical fraction as a function of 
density with their semi-analytic prescription, the mod- 
eled SO fraction does not agree with observations. Using 
a simple prescription for the relationship between the 
host dark matter halo mass accretion history and galax- 
ies lying in the halos, Maulbetsch et al. (2007) are able 
to qualitatively, but not quantitatively, reproduce the 
morphology-density and specific star formation rate en- 
vironmental correlation at ^; = 0. Additionally, the pre- 
dicted morphology-density relation predicted at z = 1 
differs dramatically from what is observed at that red- 
shift. Therefore, while it is likely that the merger history 
influences the correlations among environment, star for- 
mation history, and morphology, it is not the only driver. 

In particular, a number of baryonic processes are in- 
voked to explain these correlations. Each process is asso- 
ciated with a characteristic timescale and halo mass, pro- 
ducing different effects on the photometric and structural 
properties of either central or satellite galaxies. Central 
galaxies continue to form stars as long as they can accrete 
significant amounts of cold gas. Cosmological simula- 
tions suggest that gas can cool in halos of M < lO^^M© 



up to the current epoch, but not in more massive ha- 
los since z 2 (Keres et al. 2005; Croton et al. 2006; 
Dekel & Birnboim 2006). This trend translates directly 
into the color distributions of the associated galaxies, but 
there is no prediction for how the gas accretion history 
specifically affects morphology evolution. Satellite galax- 
ies falling into more massive halos are subject to a num- 
ber of processes, such as ram pressure stripping (Gunn 
& Gott 1972; Abadi et al. 1999; Hester 2006), "strangu- 
lation" or "starvation" (Larson et al. 1980; Bekki et al. 
2002), and harassment due to high-speed encounters with 
other galaxies (Moore et al. 1996, 1999). The relative 
importance of these mechanisms for altering morphology 
depends on both properties of the satellite galaxies and 
the larger host halo. 

Disentangling these processes is extremely difficult, 
since galaxies can experience most of them diiring some 
point in their evolution. In order to determine which pro- 
cesses drive the morphology-density, color-density, and 
star formation-density relations, it is important to ob- 
serve galaxies in a range of environments at a variety of 
epochs, in which different processes may dominate. In 
particular, since many theories on galaxy evolution fo- 
cus on processes that are most likely to be effective in 
cluster environments, where the environmental contrast 
is maximized, it is important to observe clusters in differ- 
ent stages of formation, even before the structures have 
virialized. 

Several high-redshift galaxy overdensities have been 
observed that are thought to virialize into massive clus- 
ters by the current epoch. Commonly called "protoclus- 
ters," some of these objects have been found serendipi- 
tously (Pascarelle et al.' 1996; Steidel et al. 1998, 2005). 
Others were discovered through narrow-band Lya and 
Ha surveys targeting the regions about high redshift ra- 
dio galaxies (HzRGs, Kurk et al. 2002; Overzier et al. 
2006) because HzRGs lie in highly biased environments 
(De Brcuck et al. 2000). One of the best-established (and 
most thoroughly observed) protoclusters is at z = 2.300 
in the Q1700 field centered on the HS1700-h643 quasar, 
and was discovered in a survey of ^ = 2.4 ± 0.4 star- 
forming galaxies in that field (Steidel et al. 2004, 2005). 
Spectroscopy of > 100 photometric z > 1.5 candidates in 
the Q1700 field allows this structure to be well-mapped 
in projected and redshift-space, and leads to a good es- 
timate of the mass overdensity. This stands in contrast 
to work on the other protoclusters, for which the esti- 
mated mass overdensity is based on projected number 
densities of Lyo: emitters, which do not fully trace the 
underlying population of high redshift galaxies, and for 
which the number density in average environments is not 
well quantified. Furthermore, protocluster membership 
for galaxies that are not Lya emitters is established us- 
ing photometric redshifts, which have large (5(z) > 0.5 
errors (Overzier et al. 2006). The mass overdensity of 
the Q1700 protocluster indicates that it should virialize 
into a massive cluster of ~ 1.4 x IO^^Mq by the present 
(Steidel et al. 2005). 

Stellar population synthesis modeling of the deep 
multiband photometry and spectroscopy of the star- 
forming galaxies in the Q1700 field indicate that galax- 
ies in the protocluster are, on average, significantly more 
massive and older than galaxies in the field. Complemen- 
tary analyses of the spatial clustering of z > 2 galaxies as 
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a function of rest-frame optical color suggest that redder 
galaxies are more strongly clustered (Adelberger et al. 
2005; Quadri et al. 2007), which is consistent with the re- 
sults from the Q1700 protocluster. Therefore, it appears 
that star-formation history is already a strong function 
of environment at high redshift. 

In this paper, we explore connections between mor- 
phology and environment, and morphology and star for- 
mation history at high redshift. We examine the rest-UV 
morphologies of optically-selected, spectroscopically con- 
firmed star- forming galaxies in the Q1700 field, perform- 
ing the first rigorous analysis of morphology as a func- 
tion of environment at z > 1.3. In addition, given the 
large multiband data set available for this field, we are 
able to identify a complementary population of high red- 
shift galaxies selected by their near-infrared colors, and 
compare their properties with those of the UV-selected 
galaxies. The structure of the paper is as follows: in 
Section 2, we summarize previous observations of the 
Q1700 field and present our new Hubble Space Tele- 
scope Advanced Camera for Surveys {HST / ACS) ob- 
servations, from which we derive morphologies for our 
galaxy sample. In Section 3, we describe the set of pa- 
rameters we use to quantify morphology, and present 
tests of the robustness of these parameters. We discuss 
the results of the analysis on the environmental depen- 
dence on morphology for star- forming galaxies in Section 
4. We present an analysis of morphology as a function 
of star formation history and the physical properties of 
galaxies in Section 5. In Section 6, we compare the dis- 
tributions of morphological parameters of optically and 
near- infrared selected galaxies. We discuss the results of 
our analyses and place them in the broader picture of 
structure formation and galaxy evolution in Section 7. 

For the rest of this paper, we assume a flat cosmology 
with Qrn = 0.3, = 0.7, and h = 0.7. 

2. OBSERVATIONS AND SAMPLE SELECTION 
2.1. Prior Observations 

In order to optically select high redshift star-forming 
galaxies in the foreground of the HS1700-I-643 quasar, 
Steidel et al. (2004) observed a 15^3 x 15^3 field (Q1700) 
with the Prime Focus Imager on the William Herschel 
Telescope in the Un, G, and TZ bands, with supplemen- 
tary observations using the f/„ filter with the Low Reso- 
lution Imaging Spectrometer (LRIS) on the Keck I tele- 
scope. 

Spectroscopic follow-up was concentrated on the region 
nearest the HS 1700-1-643 quasar. Redshifts were primar- 
ily measured from optical spectra probing the rest-frame 
UV. These spectra were obtained using the LRIS-B in- 
strument on the Keck I telescope (Steidel et al. 2004). 
In addition, Erb et al. (2006) measured Ha redshifts for 
19 objects in the larger optical spectroscopic sample us- 
ing the NIRSPEC near-infrared spectrograph on Keck II. 
In total, 167 galaxy photometric candidates have secure 
redshifts, of which 141 have z > 1.4. The contaminants 
were stars and low-redshift galaxies. The current Q1700 
spectroscopic sample is larger than the original one pre- 
sented in Steidel et al. (2005), due to subsequent addi- 
tional observations. 

In order to perform population synthesis modeling 
on those galaxies with secure redshifts, near- and mid- 
infrared observations of the Q1700 field were performed 



and presented by Barmby et al. (2004), Shapley et al. 
(2005) and Erb et al. (2006). Data in the 3.5 - S.O^m 
wavebands were obtained with the Infrared Array Cam- 
era (IRAC) on the Spitzer Space Telescope during the 
In-Orbit Checkout (IOC). Data reductions and photom- 
etry are described in Barmby et al. (2004) and Shapley 
et al. (2005). Using the Wide Field Imaging Camera 
(WIRC) on the 5m Palomar Hale telescope, a 8'.5 x 8'.5 
section of the Q1700 field was imaged in the J and Kg 
(Erb et al. 2006; Shapley et al. 2005) bands. These data 
also allow us to identify a sample of distant red galax- 
ies ("DRGs") satisfying the J — Kg > 2.3 criterion of 
Franx et al. (2003). This selection criterion is tuned to 
select galaxies with significant 4000A/Balmer breaks at 
z « 2.5, and hence, selects galaxies with either old or 
heavily obscured stellar populations (Forster Schreiber 
et al. 2004; van Dokkum et al. 2004; Reddy et al. 2005). 

2.2. RST/ACS Observations 

We imaged the Q1700 field with the HST /ACS us- 
ing the F814W filter. For a protocluster mean redshift 
of 2; = 2.300, the pivot wavelength of the filter corre- 
sponds to a rest-frame wavelength of A ^ 2500A, and the 
FWHM of the point-spread fimction of 0'.'12 corresponds 
to « 1 kpc at z = 2.3. In order to include as many of 
the spectroscopically confirmed galaxies as possible, we 
covered the region nearest the HS1700-I-643 quasar with 
four pointings (Figure 1). Each pointing was imaged over 
the course of five orbits of nearly equal exposure time, 
for a total exposure time of 12520 seconds. This corre- 
sponds to a sensitivity of 29.0 AB magnitude for a la 
surface brightness fluctuation in a 1 □" aperture, and 
a 28.4 AB magnitude depth for a lOcr point source in 
a O'.'l radius circular aperture. We used the MultiDriz- 
zle script (Kockemoer et al. 2002) to clean, sky-subtract, 
and drizzle the flat-fielded data products from the ACS 
CALACS software pipeline. In order to align the ACS 
image to the ground-based coordinate system, we used 
« 100 bright objects on each pointing to create a map 
between the ACS and ground-based coordinate systems. 
We then corrected the ACS image world coordinate sys- 
tem accordingly. 

2.3. Sample Selection 

High redshift photometric candidate galaxies were se- 
lected using the BX/MD/C/D/M photometric color se- 
lection schemes of Steidel et al. (2003) and Adelberger 
et al. (2004). A total of 1472 "BX", 238 "MD", 81 "C", 
74 "D", and 45 "M" candidate galaxies were identified 
in the optical ground-based observations of the Q1700 
field, for a total of 1910 optically selected galaxy candi- 
dates (Steidel et al. 2004, 2005). Of primary interest to 
us for our environmental study of the z = 2.300 galaxy 
redshift spike ("protocluster") are the BX/MD objects, 
for which the mean redshifts are z = 2.20 and z = 2.79 
respectively (Steidel et al. 2003; Adelberger et al. 2004), 
in contrast to the mean redshift z w 3 for the combined 
C/D/M sample. As discussed in Section 2.1, 141 of these 
galaxies have secure redshifts, of which 25 are in the over- 
density. In the near-infrared field, we identified 75 DRGs, 
11 of which are also classified as BX or MD objects. Of 
the DRGs that are also classified as BX or MD objects, 
7 have spectroscopic redshifts, of which two are in the 
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Fig. 1. — Schematic of observations and samples. The image covers the full scale of the 7?.-band observation. The Ks observation is 
bounded by the red box. The HST / ACS footprint is denoted by the four black-rimmed boxes. Marked by the large solid black square is the 
HS1700-I-643 quasar. Objects that satisfy the rest-frame UV color criteria but for which no redshift information is available are marked by 
small black dots; those with redshift information are denoted by large blue dots, with triangles about the spike objects. DRGs are marked 
with red squares. Objects satisfying more than one selection criterion will have stacked markers. 



overdensity. Figure 1 illustrates the footprints of all the 
observations and all of the galaxy samples. 

A total of 22 of 25 BX/MD objects whose spectro- 
scopic redshifts place them in the protocluster fell on the 
footprint of our ACS observations. Additionally, 72 op- 
tically selected objects with redshifts z > 1.7 but not 
within the protocluster were observed by ACS, as were 
234 BX/MD/ C/D/M objects without spectroscopic red- 
shifts (of which 223 are not obviously stars or large fore- 
ground galaxies) and 43 DRGs. Figures 2-5 show im- 
ages of the spectroscopically confirmed optically selected 
galaxies, and all the DRGs. All seven of the DRGs that 
also satisfy optical selection criteria and have spectro- 
scopic redshifts are on the ACS pointings, as are three of 
the four DRGs that are also BX or MD objects that do 
not have secure redshifts. A summary of the ACS galaxy 
samples can be found in Table 1. 

3. MORPHOLOGICAL PARAMETERS 



3.1. Parameter Choice 

Galaxies at z ~ 2 are not classified easily by the Hub- 
ble sequence. As demonstrated in Figures 2—5 and as 
noted in previous studies of the rest-frame UV and opti- 
cal morphologies of z > 2 galaxies (e.g., Dickinson 2000; 
Elmegreen et al. 2004b; Lotz et al. 2006; Law et al. 2007), 
high redshift galaxies tend to be clumpy and highly irreg- 
ular, with no strong evidence for spiral or elliptical types. 
Therefore, there have been several attempts to define a 
set of quantitative measures of morphology. Some (e.g., 
Ravindranath et al. 2006; Zirm et al. 2006) have used 
circularly symmetric Sersic fits to quantify the surface 
brightness distribution. Others have used the CAS (con- 
centration C, asymmetry A, dumpiness S) system (Kent 
1985; Schade et al. 1995; Bershady et al. 2000; Conselice 
2003). However, both the Sersic fits and the C and A 
parameters of CAS are defined assuming that elliptical 
isophotes provide a good description of the objects in 
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Fig. 2. — HST /ACS images of galaxies spcctroscopically confirmed to be in the protocluster, sorted in order of increasing gini. Each 
thumbnail is 4" X 4" in size. Images arc centered at the centroid of the 7^-band isophotes. Tiie elongated morphology of BX 913 may be 
due to do lensing by the foreground X-ray cluster RXJ1701. 3+6414 at 2 = 0.453 (MuUis et al. 2003). 

to estimate galaxy morphologies. These coefRcients have 
the advantage of not enforcing symmetries on the galaxy 
light distribution. The parameters we use are: size, gini, 
and multiplicity. 



TABLE 1 
Galaxy Samples 



Sample 



Number[SEDs]' 



BX/MD spike 


22 


15 




BX/MD non-spike 


63 


45 




total BX/MD w/redshift 


85 


60 





BX/MD no redshift 



197 



total BX/MD 
DRG 

LBG (C/D/M) with redshift 
LBG (C/D/M) no redshift 
total LBG (C/D/M) 



282 

43'' 

9[2] 
26 
35 



'^SED fits were obtained for objects with near- 
and mid-IR imaging and spectroscopic redshift 
information. These are the objects used in the 
comparison of physical parameters and morphol- 
ogy- 

''Some DRGs are also selected optically (10 to- 
tal) and have spectroscopic redshifts (7 of 10); 
however, redshifts are not available for most of 
the sample. 



3.1.1. Size 

The simplest way to quantify the size of a galaxy at 
high redshift is by its projected (angular or physical) 
surface area. In particular, since the angular diameter 
distance in our cosmology is < 10% smaller at the upper 
end of the BX/MD spectroscopic sample redshift distri- 
bution than at the lower end, we use the angular pro- 
jected surface area as the measure of galaxy size. This is 
also a pragmatic choice given that we do not have spec- 
troscopic redshifts for some of our galaxy samples. For 
simplicity, we define size 



I = N, 



pixi 



(1) 

where Npi^ is the number of pixels associated with the 
galaxy. To convert to angular units, one must multiply / 
by the angular size of an ACS pixel ((0'.'05)^). A typical 
size of 250 ACS pixels corresponds to 0.625 □", or 43 
kpc^ at ^ = 2.3. 



question. These parameters arc also extremely sensi- 
tive to the precise location of the "center" of the galaxy. 
Given the highly irregular nature of high redshift galaxy 
morphologies, it is unclear if elliptical, centroid-sensitive 
fits describe the galaxies in a robust and meaningful way. 
Instead, we adopt a set of non-parametric coefficients 



3.1.2. Gini 

An important measure of galaxy morphology is the dis- 
tribution of light emission, in the sense that one would 
like to know if the galaxy surface brightness is roughly 
constant, or if there are large variations in surface bright- 
ness. For galaxies with a high degree of symmetry (for 
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Fic;. 3. HST / ACS images of galaxies spectroscopically confirmed not to be in the protocluster. As in Figure 2, galaxies are sorted in 
order of increasing giiii, sijamiiiig the range of 0.245 to 0.449. Images are centered on the centroids of the 7^-band isophotes. See Figure 4 
for ginis in range 0.449 to 0.790. 



example, galaxies that are easily classified by the Hubble 
sequence), measures such as the Sersic index or the con- 
centration parameter C (proportional to the logarithm 
of the ratio of the radius containing 80% of the light in 
the galaxy to the radius containing the inner 20%) in- 
dicate the strength of variations in surface brightness. 
However, since high rcdshift galaxies are irregular, we 
instead adopt the gini parameter 



G = 



Y^{2i-Npi,-l)Xi, (2) 



where X is the mean sky-subtracted flux per pixel, Xi 
is the sky-subtracted flux in pixel i, and the pixels have 
been sorted in order of increasing flux. This parameter, 



originally used in economics to describe the distribution 

of wealth in a population, was first used in an astronom- 
ical context by Abraham et al. (2003), and its usefulness 
in the context of galaxy morphology has been examined 
by Lotz et al. (2004, 2006) and Law et al. (2007). 

Gini is defined to have a value G = if the flux is uni- 
formly distributed among pixels, and G = 1 if one one 
pixel hoards all the flux in the galaxy. Thumbnails of the 
spectroscopic BX/MD and the DRG samples are sorted 
in order of increasing gini (Figures 2-5). Note that ob- 
jects with the same gini may span the range of diffuse 
objects with multiple small knots to those with a single 
bright central region. Therefore, unlike the Sersic index 
or C, gini can accommodate a wide range of spatial light 
distributions. 
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Fig. 4. — Continuation of Figure 3. HST/ ACS images of galaxies spectroscopically confirmed to NOT be in the redshift spike. Galaxies 
are sorted in order of increasing gini, spanning the range 0.449 < G < 0.790. See Figure 3 for galaxies with ginis in the range 0.245 to 
0.449. 



3.1.3. Multiplicity 

In addition to describing how light is distributed 
among pixels, one would like to characterize the shape 
or dumpiness of the galaxy light distribution. We use 
the multiplicity ('5) parameter introduced by Law ct al. 
(2007) as a measure of this property. Multiplicity is anal- 
ogous to the total potential energy of a gravitational sys- 
tem. The total "potential energy" of the flux distribution 
can be described as 



were in the most compact configuration (i.e., with the 

brightest pixel at the center, and the next brightest pix- 
els forming a ring around that pixel, and so forth), this 
configuration would have the greatest "total potential 
energy" : 



V'. 



compact 



E E 

1=1 j=ij>^i 



(4) 



E E 

j=i.j5^i 



XiXj 



(3) 



where r'^j is the distance between the ith and jth particles 
in the new configuration. The multiplicity is then defined 
to be 

^compact 



* = 100 logio 



where rij is the distance (in pixels) between the ith and 
jth pixels. If one were to reorder the pixels so that they 



(5) 



The range of ^ depends somewhat on the method used 
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Fig. 5. — HST / ACS images of objects satisfying the J — Kg > 2.3 DRG color criterion. Objects are sorted in order of increasing gini. 
Images are centered at the centroid of the ifg-band isophotes. 
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Fig. 6.— Sample galaxies (MD 109, BX 939, BX 575, and BX 
592) demonstrate the range and interpretation of multiplicity. 

to assign pixels to galaxies (see Section 3.2). Generally, 
objects with < 2 arc bulgc-likc and objects with < 5 
tend to be dominated by one main clump. Objects with 
high multiplicity > 15) tend to have many almost 
equally bright sources that are widely separated. Galax- 
ies with intermediate multiplicity can be highly elon- 
gated, or have several clumps that arc close together. In 
general, galaxies that consist of several clumps of uneven 
luminosity will have lower multiplicity than a galaxy for 
which all the clumps are of approximately equal lumi- 
nosity. Figure 6 demonstrates galaxies that are typical 
of their multiplicities. 

3.2. Pixel Selection 

A major problem in characterizing the morphologies 
of high redshift galaxies is determining which pixels arc 
associated with the galaxy, and how the pixel selection 
affects the quantitative measures of morphology. Lotz 
et al. (2004) select pixels by first determining the ellip- 
ticity and Pctrosian radius of a galaxy, using elliptical 
apertures. They then convolve their data within the Pct- 
rosian ellipse with a gaussian, and select all pixels within 
the Petrosian ellipse that meet a surface brightness cut- 
off. As discussed in Law et al. (2007), the caveats associ- 
ated with this approach include the fact that it assumes 
symmetry for a population of highly irregular objects, 
and it tends to associate many sky pixels with the galaxy. 

Law et al. (2007) use a method that does not depend on 
the radial symmetry of the object. An initial estimate of 
the galaxy centroid is obtained from ground-based imag- 
ing. The galaxy centroid is then recalculated, based on 
the light distribution in the ACS images. Next, an aper- 
ture of I'.'S radius is used to search for pixels that are 
associated with the galaxy. Pixels are assigned with a 
galaxy if they meet a redshift-dependent surface bright- 
ness threshold. The redshift-dependent surface bright- 
ness threshold scales with cosmological dimming, and 
its purpose is to prevent bias in morphological measure- 
ments. This threshold is defined as n{z)a, where a is the 
standard deviation in the sky background. Law et al. 
(2007) set the threshold for the highest redshift galaxy 
to n{zhigh) = 3, where Zhigh = 3.4. Then, since surface 
brightness dims as (1 + z)~^ for a fixed passband, they 
set the surface brightness threshold 



n{z) = 3 



1 



Zhigh 



1 3 



1 + Z 



(6) 



The advantage of this surface brightness threshold 
method is that the geometry of the object does not mat- 
ter. This disadvantage of this method, for our purposes, 
is that that it requires redshift information, which is 
missing for many of our galaxies, in particular, the ma- 
jority of the DRGs. 



For the current analysis, we adopt a pixel selection 
method that is redshift independent and includes all the 
pixels associated to the galaxy without including sky pix- 
els. We use a redshift independent selection method 
for two reasons. First, given the poor constraints on 
photometric redshifts, it does not make sense to use a 
redshift-dependent threshold based on photometric red- 
shifts for our samples without spectroscopic information. 
Secondly, the redshifts of our spectroscopic sample span 
a smaller range than those of Law et al. (2007) and Lotz 
et al. (2006). Lotz et al. examine two samples with mean 
redshifts z ^ 1.5 and z ^ 4, while Law et al. break up 
their 1.8 < z < 3.4 sample into two samples with mean 
redshifts of z ^ 2 and z ^ 3. More than 75% of the galax- 
ies in our BX/MD sample fall into the lower-redshift bin 
from Law et al. (2007). Furthermore, though they span 
a larger redshift range than confirmed protochister mem- 
bers, galaxies in the Q1700 field that are not contained 
within the redshift spike have a similar mean redshift, 
z = 2.34. The similarity in the mean redshifts is due to 
the fact that z = 2.3 lies near the peak of the BX/MD 
redshift selection function (Stcidcl et al. 2004). There- 
fore, it is less crucial in the current study to account for 
cosmological dimming in the surface brightness thresh- 
old. 

The first step in selecting pixels associated with a par- 
ticular galaxy is to define the region in which fiux could 
reasonably be associated with the galaxy. Like Law et al. 
(2007), we begin by recalculating the galaxy centroid in 
order to center the pixel selection region. The centroid is 
calculated based on the first-order moment of light on the 
ACS images, usually calculated within a 30-pixel (1'.'5) 
radius circle centered on the TZ centroid. We use a differ- 
ent aperture within which to calculate the centroid if the 
multiband ground-based images indicate that there is an- 
other object (either a large foreground galaxy or galaxy 
with different colors than the object in question) within 
the default 1'.'5 aperture, or if the 7?,-band isophote has 
a linear extent > 3". In those cases, the size of the 
centroid-finding aperture is tailored to avoid contami- 
nants and completely surround the 7^-band detection. 
In some cases, a single ground-based detection will en- 
compass several clumps of light with separations of ^ 1" 
(i.e., roughly the size of the ground-based seeing), which 
corresponds to a projected physical separation of 8 kpc 
at z = 2.3. It is highly unlikely that such small separa- 
tions are the result of a projection effect, so we consider 
all clumps that fall within the ground-based detection 
isophote to be associated with each other. Therefore, we 
calculate the centroid for all patches of light that lie in 
the ground-based detection. 

Next, we define an aperture about the newly calculated 
centroid within which we search for pixels associated with 
the galaxy. As in Law et al. (2007), the default aperture 
size is a 1'.'5 radius circle. However, as in our search 
for the galaxy centroid, we must customize the aperture 
size in some cases. In those cases, we use the ground- 
based 7^-band images as a guide, since the seeing-limited 
7?,-band detections provide an approximate upper limit 
to the extent of the galaxy. If the 7?.-band detection 
extends beyond I'.'S (~ 12 kpc in projection) from the 
center of the object, we use a circular aperture on the 
ACS image that encompasses the 7?.-band detection. If 
there are foreground galaxies or objects that have signifi- 



10 



Morphologies of Protocluster Galaxies 




Fig. 7. Pixel selection for BX 759 (z = 2.418), which would be 
classified as "nebulous" in Law et al. (2006). Lejt: BX 759 with- 
out pixel selections marked. Image size is 2" x 2". Right: White 
contours mark pixel selection according to the redshift-dependent 
pixel selection scheme of Law et al. (2007; Scheme B). This cor- 
responds to selecting all pixels > Sit above the sky at x = 2.418 
such that galaxies at z/tjg/t = 3.4 are selected with a Za surface 
brightness threshold. Black lines delineate the regions which are 
associated with the galaxy using the fiducial pixel selection scheme. 

cantly different colors than the object in question that lie 
within the default aperture, which occurs for ~ 15% of 

the spectroscopic sample, wc use circular apertures that 
are small enough to avoid contaminants. It was only 
necessary to use a non-circular aperture twice (for the 
photometric BX/MD sample) in order to avoid contami- 
nants. In these cases, the elliptical apertures were aligned 
to the long-axis of the galaxy and made large enough to 
encompass the whole galaxy, but small enough to avoid 
contaminants. The size of those elliptical apertures was 
also adjusted so that the total area within the aperture 
was not signifcantly different from the area of the default 
aperture. 

As in Law et al. (2007), we perform local sky- 
subtraction for each object to account for any imperfec- 
tions in the MultiDrizzle sky subtraction. We calculate 
the sky background in an annulus about the detection 
aperture, typically 6 □" (or ^2400 pixels) in total sur- 
face area. We use a sigma-clipping algorithm to exclude 
any pixels that vary by more than Str about the sky level 
from the calculation of the sky background. 

In order to determine which pixels in the detection 
aperture are associated with the galaxy, we convolve the 
data with a circularly symmetric gaussian with a one- 
pixel standard deviation, truncating the gaussian four 
standard deviations from the center, using the IRAF tool 
"gauss." This allows us to exclude sky pixels from our 
segmentation map, since smoothing reduces the signif- 
icance of isolated bright pixels. Wc select pixels to be 
associated with the galaxy if the smoothed flux is more 
than 2fT above the sky, with cr calculated from the raw 
data. We then use the raw sky-subtracted flux in those 
pixels for our morphological analysis. In the rest of the 
paper, we refer to this pixel selection scheme as the "fidu- 
cial" scheme. Figure 7 shows pixels selected with both 
our fiducial method and the method of Law et al. (2007). 

In the interest of determining the robustness of mor- 
phological parameters derived using our pixel selection 
scheme, we compare the values of size, gini, multiplicity, 
and the total flux using our pixel selection with results 
from various other pixel selection methods (Figure 8). 
Points on Figure 8 correspond to all high redshift galaxies 
with spectroscopic redshifts, including C/D/M galaxies, 
~ 3. We consider three reasonable schemes for select- 
ing pixels. In Scheme A, all pixels with fluxes satisfying a 



5(T threshold on the unsmoothed imaged within a 1'.'5 ra- 
dius of the ccntroid are associated with the galaxy. This 
allows us to evaluate the effects of smoothing for con- 
stant surface brightness selection schemes. In Scheme B, 
we use the redshift-dependent pixel selection scheme of 
Law et al. (2007) to associate pixels with galaxies. This 
allows us to examine the effects of a redshift-dependent 
surface brightness threshold. In Scheme C, we use the 
fiducial scheme with the one modification that we fix the 
detection aperture to 1'.'5, so that we can isolate the ef- 
fects of aperture size on the morphological parameters. 
In Figure 8, the dotted line in each panel is the line of 
equality {y = x), and the solid line is the least-squares 
fit to the data. 

Since we are interested in the differences of the mor- 
phologies of different samples (spectroscopically con- 
firmed protocluster and "field" BX/MDs, BX/MDs with- 
out redshifts, and DRGs). it is less important that the 
quantitative measures of morphology be absolutely cali- 
brated as a function of pixel selection than it is that they 
be robust in a differential way. In other words, we would 
like a strong correlation to exist between morphological 
parameters calculated with two different pixel selection 
schemes, but it is not necessary that the parameters take 
the same value. We see that the sizes measured with dif- 
ferent pixel selection schemes are well-correlated (see Fig- 
ure 8), as are the ginis and fluxes. This is especially true 
for Scheme C, since it differs from the fiducial scheme 
only in the size of the detection aperture. 

Multiplicity shows much more scatter than the other 
parameters. However, most of the scatter comes from 
only two effects. Points in Figure 8 outlined with blue 
diamonds represent galaxies for which the 1'.'5 detection 
aperture includes contamination from other galaxies that 
are nearby in projection. Because multiplicity is sensitive 
to the spatial distribution of light, the effect of the con- 
taminants is to systematically increase the multiplicity. 
Points outlined with red squares indicate galaxies that 
have G < 0.15 as calculated in Scheme B (i.e., in the fash- 
ion of Law et al. 2007) but have high multiplicity. Law 
et al. (2007) classify galaxies with these parameters as 
"nebulous," and comment that multiplicity is ill-defined 
for these galaxies. Since they tend to have low surface 
brightness, they are quite sensitive to the pixel selection 
method. We have designed our fiducial scheme to cap- 
ture as much information from these objects as possible 
(see Figure 7 for an example of which pixels the fiducial 
scheme and Scheme B select for a "nebulous" galaxy). 

Therefore, we can conclude that, for the most part, our 
set of morphological parameters is robust in a differential 
sense to pixel selection. We use our fiducial pixel selec- 
tion method throughout the paper, but will occasionally 
check results from the fiducial scheme against those from 
the other schemes. 

Even for a fixed pixel selection technique, there are un- 
certainties in morphological parameters due to the finite 
signal-to-noise (S/N) of the object detections. We assess 
such uncertainties by running simulations of our morpho- 
logical analysis procedure with the fiducial pixel selection 
technique. Using pixel flux distributions drawn from our 
actual sample of objects, spanning a range of size, gini, 
and multiplicity, we add copies of each object light distri- 
bution to 100 random locations in the ACS image, and re- 
cover the morphological parameters for simulated objects 
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Fig. 8. — Morphological parameters as a function of pixel selection. Blue diamonds mark galaxies for which the 30 pixel aperture captures 
flux from nearby objects. Red squares mark galaxies that classified as "nebulous" by Law et al. (2007, G < 0.15 in their scheme) for which 
multiplicity is ill-defined. 



in relatively clean locations (i.e. objects whose fluxes are 
not significantly contaminated by nearby bright objects). 
This procedure indicates that the typical size uncertain- 
ties are 5 — 10%, the typical uncertainties on gini values 
are < 5% (i.e. a(G) = 0.01 - 0.02), and the typical mul- 
tiplicity uncertainties are 5 — 10%. For the faint DRGs 
in our sample, the uncertainties are larger, at the level 
of ~ 15% for size, 10% for gini (i.e. a[G) = 0.02), and 
30% for multiplicity. 

4. SPIKE VS. NON-SPIKE OBJECTS: MORPHOLOGY AS A 
FUNCTION OF ENVIRONMENT 

Our primary interest in determining if there is a de- 
pendence of morphology on environment comes from the 
initial analysis of the Q1700 spike, which showed that cer- 
tain properties of the spike galaxies were quite different 
from those not in the spike (Steidel et al. 2005). Empir- 
ically, galaxies within the spike have redder TZ — Kg col- 
ors than galaxies outside of the spike. Stellar population 
synthesis modeling of the broadband colors of the spec- 
troscopically confirmed BX/MD objects showed that, re- 
gardless of whether models were fit with a constant or 
exponentially declining star formation rate, galaxies con- 
tained in the spike are significantly older and more mas- 
sive than those not in the spike. These facts suggest that 
star formation history is already a function of environ- 
ment at z = 2.3. 

In this section, we will consider morphology as a func- 
tion of environment; in the next section, we consider mor- 
phology as a function of physical properties (mass, age, 
etc.) and star formation history. In this way, we can de- 
termine if environment, star formation history, and mor- 
phology are as tightly related at high redshift as they 
are in the current epoch. As a cautionary note, we em- 



phasize that our data probe rest-frame UV, not optical 
wavelengths, and so we are probing morphologies of the 
young star-forming regions of each galaxy instead of the 
stellar mass as a whole. However, there is some evidence 
(Dickinson 2000; Zirm et al. 2006) that the rest-frame 
UV morphology is a reasonable tracer of optical mor- 
phology (but see Toft et al. 2005). 

For the analysis of the morphology distributions as 
a function of environment, we must first define what 
is meant by "environment." In large surveys, such as 
SDSS and the DEEP2 galaxy redshift survey (Davis et al. 
2003), "environment" refers to the number density of 
galaxies, as calculated on Mpc scales. In our context, 
environment refers not to the local number density of 
galaxies, but to the binary designation of whether or not 
a galaxy belongs to the protocluster. We define environ- 
ment in this way due to our small sample size. Therefore, 
we separate our spectroscopic BX/MD sample into two 
subsamples. The "spike" sample consists of 22 BX/MD 
galaxies with spectroscopic redshifts of 2; = 2.300±0.015. 
As shown in Figure 1 of Steidel et al. (2005), the num- 
ber of galaxies in this small redshift interval exceeds the 
average expected number according to the BX/MD red- 
shift selection function, by a factor of ~ 7. Typical 
redshift uncertainties of az ~ 0.002 (Adelberger et al. 
2005) do not significantly bias the estimate of the spike 
overdensity. The "non-spike" sample consists of 63 spec- 
troscopically confirmed BX/MD galaxies with z = 2.34. 
Since we are interested in comparing the morphologies 
of galaxies at the same epoch but in different environ- 
ments, we do not include the additional C/D/M galaxies 
because their high mean redshift z — 2.94. The spec- 
troscopic sample of non-spike BX/MDs serves as a good 
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control sample due to the fact that its mean redshift is 
almost the same as that of the spike. We calciilate the 
size, gini, multiplicity, and total flux of each galaxy using 
the fiducial pixel selection method of Section 3. 

To evaluate the differences between the spike and 
non-spike samples in the distributions of morphologi- 
cal parameters (Figure 9), we perform a one-dimensional 
Kolmogorov-Smirnov (KS) test for each parameter. As 
demonstrated in Figure 9 and Table 2, there is no statis- 
tically significant difference between the spike and non- 
spike distributions for each parameter. In order to test 
the robustness of these findings, we perform KS tests for 
a variety of pixel selection methods. Also listed in Table 
2 are the KS probabilities that the spike and non-spike 
morphologies are drawn from the same distribution as 
calculated for Schemes A, B, and C. For no pixel selec- 
tion method are the spike morphological parameter dis- 
tributions significantly different from the non-spike dis- 
tributions. 

We also perform an analysis to determine if there are 
different proportions of "bulge-dominated" galaxies in 
the two samples, since, in the local universe, bulge- 
dominated systems (ellipticals and SO's) are preferen- 
tially found in dense environments. We define a galaxy 
to be bulge- like if has a multiplicity < 2 and G > 0.5, 
which means that it is a single compact system with a 
sharply falling flux distribution function. While objects 
with < 2 are all spatially compact, they can have 
flux distribution functions characteristic of spheroids, ex- 
ponential disks, or something in between. The cut in 
gini separates galaxies with steep flux distribution func- 
tions (which would be characteristic of bulge-dominated 
systems) with those that have flux more uniformly dis- 
tributed across pixels (typical of disk- like galaxies). Sep- 
arating the bulge population using cuts in two of our mor- 
phological parameters is analogous to the classification of 
bulges in Lotz et al. (2006) using gini and the second mo- 
ment of the light distribution. This two-dimensional clas- 
sification of bulges is a cruder metric than, for example, 
the Sersic index, and encodes no dynamical information. 
However, it is more useful if the bulk of the galaxy sam- 
ple has irregular morphology, and does separate galaxies 
that appear to be compact and roughly axially symmetric 
in projection from the more irregular galaxies. Using this 
metric, we find that 3 of the 22 (14%) spike galaxies are 
bulge-like, as are 11 of the 63 (17%) non-spike objects. 
Given the small number of galaxies in each sample, the 
difference in bulge fraction is not statistically significant. 

In order to explain the lack of differences in morphol- 
ogy between the spike and non-spike samples, there are 
four main points of interest. The first is that our sample 
size is fairly small. One can detect large differences in dis- 
tributions relatively easily with small samples, for exam- 
ple, the large differences in the mean stellar mass of the 
spike and non-spike galaxy samples. However, it would 
be difficult to ascertain subtle differences with our small 
samples using the KS test, especially if the differences 
were mainly in the tail of the distributions, for which the 
KS test is not the optimal statistic. The second point is 
that, with the F814W filter, we are probing a rest-frame 
wavelength that is in the UV. If we are interested in mor- 
phology as a proxy for stellar mass distribution, it would 
be better to probe the galaxies at a rest-frame wave- 
length that is less sensitive to the brightest (youngest) 



TABLE 2 

KS Test Comparison of Morphology and Environment 



Pixel Selection Scheme 


Size 


Gini 


Multiplicity 


Flux 


Fiducial 


0.78* 


0.61'' 


0.18'' 


0.69* 


Scheme A 


0.91 


0.73 


0.59 


0.88 


Scheme B 


0.27 


0.33 


0.83 


0.44 


Scheme C 


0.68 


0.58 


0.71 


0.54 



''The KS probability that the spike and non-spike morphologies 
are drawn from the same distribution. 



stars and dust extinction. While several groups (e.g., 
Dickinson 2000; Zirm et al. 2006) claim that differences in 
rest-frame UV and optical morphologies are small, these 
conclusions are based on very small samples imaged with 
the NIC3 camera on HST, which has a coarser resolution 
( 0'.'26) than that of ACS. Thirdly, we have analyzed 
the morphologies of galaxies that were photometrically 
selected to be star-forming and relatively unobscured. 
Hence, we are missing populations of old or obscured 
galaxies from our analysis; one might worry that we are 
doing the equivalent of measuring the morphologies of 
only spiral and irregular star-forming galaxies in the local 
universe. Therefore, according to these arguments, there 
could still be real differences in the spike and non-spike 
morphology distributions even though we do not detect 
significant evidence of differentiation with our analysis. 

The fourth point is that, since the galaxies must be 
young (as indicated by their high redshift and stellar 
ages, typically ~ 100 Myr — 1 Gyr ), they may not have 
had time to settle into dynamical equilibrium. Stud- 
ies of galaxies forming by monolithic collapse (Katz & 
Gunn 1991; Katz 1992; Immeli et al. 2004) or by mergers 
(Robertson et al. 2004) demonstrate that the timescales 
for galaxies to reach a recognizable disk or bulge mor- 
phology are of order Gyrs. Therefore, morphology may 
be decoupled from star formation history at high redshift. 
We explore this connection in the following section. 

5. MORPHOLOGY AND PHYSICAL PROPERTIES OF 
GALAXIES 

In order to examine morphology in the context 
of the star formation history and physical properties 
of galaxies, we use properties derived from the stel- 
lar population synthesis modeling of the broadband 
UnGTZJKs+Spitzer/lRAC photometry performed by 
Shapley et al. (2005) and Erb et al. (2006). These pa- 
rameters include stellar mass, age, extinction, and star 
formation rate. In the case of Erb et al. (2006), Kg 
band magnitudes are corrected for the Ho; emission line. 
Stellar population parameters were only estimated for 
galaxies with both spectroscopic redshifts and near-IR 
photometry. The criterion on the redshifts was imposed 
because simultaneous fitting of photometric redshifts and 
star formation history leads to large uncertainties in both 
the photometric redshifts (typically \zspec — Zphot\ = 0.3) 
and stellar population model parameters (Shapley et al. 
2005; van Dokkum et al. 2006). It is necessary to have 
near-IR photometry because near-IR bands lie on the 
red side of the age-sensitive 4000A and Balmer break 
features, and hence, are crucial in separating young from 
old stellar populations. There is a total of 60 BX/MD 
objects on the ACS pointings that satisfy both observa- 
tional requirements, including 22 spike galaxies. 
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Fig. 9. — Distributions of morphological parameters (size, gini, multiplicity, flux). The morphology distributions of spike objects are 
denoted with magenta solid lines. The distributions of the non-spike galaxies are denoted with cyan hatched lines. The KS probability 
that the distributions of spike and non-spike objects are drawn from the same distribution is located in the upper right-hand corner of each 
plot. 



5.1. Methods: Stellar Population Synthesis Models and 
Statistical Tests 

The details of population synthesis modeling are de- 
scribed in Shapley et al. (2005) and Erb et al. (2006). 
Briefly, theoretical SEDs are derived assuming solar 
metallicity and a Chabrier (2003) initial mass func- 
tion (Bruzual & Chariot 2003). Starlight is attenuated 
with the Calzetti et al. (2000) dust attenuation model, 
which can be parameterized in terms of the extinction 
E{B — V). The star formation history is assumed to be 
a single burst t— model, such that the star formation rate 
SFR{tsf) = SFRoe-'^^f/'', where SFRq is the star for- 
mation rate at the beginning of the star-forming epoch, 
tsf is the time since the beginning of star formation (i.e., 
it is the age of the stellar population) , and r parameter- 
izes the length of the epoch of star formation. Constant 
star formation is equivalent to the statement that r = oo. 
The for each model is computed by taking the differ- 
ence between the broadband photometry and the model 



SED. 

There are several items to note for our analysis below. 
First, each model is parameterized by stellar age, stellar 
mass, star formation rate (SFR), extinction E{B — V), 
and T. Secondly, given the degeneracy between dust ex- 
tinction and the age of the stellar population (or more 
generally, the star formation history) , the distribution 
is fairly broad for most galaxies. Fits with a constant star 
formation rate (r = oo) are generally almost as good as 
the best fit model. This means that, in general, stellar 
age, E{B — V), SFR, and r are loosely constrained. How- 
ever, many authors (Papovich et al. 2001; Shapley et al. 
2001, 2005) have shown that the assembled stellar mass 
is more robustly constrained. Below, we will consider 
the model parameters of both the best fit stellar popula- 
tion synthesis model and the best fit model assuming a 
constant star formation rate. 

In addition to examining the relationships between star 
formation history and morphology, we also investigate 
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TABLE 3 

Correlations of morphological and physical parameters for BX/MD objects 







Size 

KS^' 




Gini 

KS 




Multiplicity 
KS N,^ 


Constant SFR 


Age 


0.059142 


1.20 


0.275269 


1.92 


0.275269 


-0.75 




Mass 


0.059142 


-1.48 


0.023213 


1.87 


0.275269 


-0.96 




E(B-V) 


0.965485 


-0.49 


0.023213 


1.19 


0.059142 


1.03 




SFR 


0.275269 


-1.53 


0.770953 


-0.90 


0.134947 


1.09 


Best Fit 


Age 


0.134947 


1.41 


0.497342 


1.74 


0.134947 


-1.08 




Mass 


0.023213 


-1.69 


0.023213 


1.66 


0.134947 


-1.39 




E(B-V) 


0.770953 


0.79 


0.134947 


1.01 


0.002571 


1.97 




SFR 


0.965485 


1.09 


0.965485 


-0.67 


0.275269 


1.20 


UV Absolute Magnitude 




0.002929 


3.73 


< 0.000001 


5.29 


0.490256 


-0.33 



''The KS probability that the morphological parameters in the lowest physical parameter bin are 
drawn from the same distribution as those in the highest physical parameter bin. 
'^Number of standard deviations from the null hypothesis for the Spearman correlation test. 
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Fig. 10. — The correlations between stellar age (derived from the best-fit stellar population model with constant star formation rate) 
and morphological parameters. Open triangles mark spike galaxies. Left: Size as a function of age. Right: Gini as a function of age. 



any connections between morphology and rest-frame UV 
luminosity. To determine the rest-frame UV luminosity, 
we use the G-band magnitude (centered at 4780 A) and 
neglect K-corrections and extinction corrections. There- 
fore, the UV luminosity is centered at wavelengths de- 
pendent upon the redshift of the object; the G band cor- 
responds to a rest-frame luminosity centered at 1450A 
for the spike mean redshift of z = 2.300. For the range 
of redshifts including 80% of the BX/MD spectroscopic 
sample {z = 1.9 — 2.8), the G band probes rest-frame 
luminosities centered at 1650 — 1260A. 

We used two different statistical measures to evaluate 
relationships between morphological and physical prop- 
erties. For the first test, we sorted galaxies by the phys- 
ical parameter in question. We then divided the galax- 
ies into three bins, and performed one-dimensional KS 
tests on the morphological distributions in the highest 
and lowest bins (as determined by the physical parame- 
ter). This test allows us to explore if, for example, the 



most massive third of galaxies have different morphol- 
ogy distributions than the least massive third of galax- 
ies. Secondly, we performed Spearman rank correlation 
tests for each pair of morphological and physical param- 
eter distributions. The results of our analysis are shown 
in Table 3. We quantify the results of the KS tests in 
terms of the probability that the morphological distribu- 
tions in the uppermost and lowermost physical parameter 
bins are drawn from the same distribution. We quantify 
the results of the Spearman test in terms of the number 
of standard deviations from the expectation value of the 
hypothesis that the data are uncorrelated. 

5.2. Morphology and Stellar Populations 

The main result from our statistical tests is that there 
are very few strong correlations between physical proper- 
ties and morphology. First, we consider stellar mass and 
age, which in Steidel et al. (2005) showed significant dif- 
ferences in distribution as a function of environment. For 
both the best fit r and constant star formation models, 
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Fig. 11. — The correlations between stellar mass (derived from the best-fit stellar population model with constant star formation rate) 
and morphological parameters. Open triangles mark spike galaxies. Left: Size as s function of mass. Right: Gini as a function of mass. 



age and mass are only weakly correlated with galaxy size 
(at the ~ 1.5ct level) and gini (at slightly less than 2a). 
In Figures 10 and 11, we show scatter plots in the age- 
size, age-gini, mass-size, and mass-gini planes to demon- 
strate how tenuous the correlations are. Moreover, our 
statistical tests demonstrate that multiplicity is not well 
correlated with either stellar age or stellar mass. The 
fact that morphology correlates only weakly with stellar 
age and stellar mass reconciles the finding that while age 
and mass are strongly correlated with environment (Stei- 
del et al. 2005), morphology is not. 

Next, we comment on the extinction and star forma- 
tion rate results. As demonstrated in Table 3, the cor- 
relation strength between morphology and E{B — V) or 
SFR are generally weak. Moreover, the deviations from 
the null hypothesis that morphology and E{B~V) or the 
SFR are uncorrelated are not always in the same direc- 
tion, depending on whether the physical quantities are 
estimated from the best fit r-model or constant star for- 
mation model. The lack of consistency is a statement 
about the inability of the stellar population synthesis 
models to strongly constrain the values of E{B — V) 
and SFR. Given the large error bars on these parame- 
ters, it is also unlikely that the KS and Spearman tests 
yield believable information on the strength of correla- 
tion between the E{B — V) and morphology, and SFR 
and morphology. It should be noted that our results are 
in contrast with those of Law et al. (2007), who find a 
significant correlation between E{B — V) and both size 
and gini. Law et al. (2007) estimated E{B — V) using 
best-fit r-models. Given the uncertainty in estimating 
the extinction with SED fitting, especially when r is al- 
lowed to vary, it is not so surprising that the correlation 
tests for these relatively small samples yield different re- 
sults. 



5.3. Morphology and UV Luminosity 



We find the strongest correlations between UV lumi- 
nosity (or absolute magnitude) and the morphological 
parameters gini and size (Figure 12, Table 3), while there 
is no evidence for a correlation between multiplicity and 
UV luminosity. We perform several tests to determine if 
the observed correlations with UV luminosity are robust. 
First, we perform Spearman rank correlation tests for a 
variety of pixel selection schemes, to ascertain that the 
correlations do not depend on the pixel selection method. 
As demonstrated in Table 4, the correlations appear to 
be robust to pixel selection method. 

Secondly, we consider that surface brightness depen- 
dent pixel selection schemes might bias the luminosity- 
size and luminosity-gini relationships. The fundamen- 
tal reason why one ought to be concerned about surface 
brightness dependent assignments to galaxies has its ori- 
gin in our definitions of gini and galaxy size. Gini is a 
parameter that quantifies the flux distribution function 
(i.e., the number of pixels of given fiux). However, given 
that our instrument does not have infinite sensitivity, we 
calculate gini with only those pixels whose surface bright- 
ness is above a certain threshold. To make an analogy 
with economics (the field which gave birth to the gini 
coefficient), applying a surface brightness threshold cut 
would be the economic equivalent of including only those 
whose income is above the poverty line in calculating the 
distribution of wealth in a population. In economic lit- 
erature, such a restriction on the distribution is called a 
"truncation bias." Since we also define size in a surface 
brightness-dependent way, as opposed to the traditional 
half-light radius, size is also subject to truncation bias. 

It is important to consider what a constant surface 
brightness threshold means when considering galaxies 
with a range of luminosities. Using an economic analogy, 
we want to consider what gini means if we were calcu- 
lating gini for two countries (e.g., the United States and 
Mexico) with different GDPs, using incomes measured 
above a threshold fixed to the poverty line of one of the 
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Fig. 12. — The correlation between UV absolute magnitude and galaxy size (left) and gini (right) for the spectroscopic BX/MD sample. 
Open triangles mark spike galaxies. 



TABLE 4 

Correlations with UV absolute magnitude 



Scheme 


Size 


Gini 


Multiplicity 


Flux 


Fiducial 


3.73=' 


5.29=* 


-0.33=^ 


7.50'' 


A 


5.40 


5.17 


-1.45 


7.17 


B 


6.39 


5.49 


-1.07 


7.27 


C 


4.05 


5.24 


-0.23 


7.26 



''Number of standard deviations from the null 
hypothesis that the given parameter is uncorre- 
lated with UV absolute magnitude, as determined 
by the Spearman rank correlation test. 
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Fig. 13. — Schematic flux distribution functions for two galaxies 
of different luminosity but with the same form of the distribution 
function. 



countries (e.g., fixing the threshold for both countries at 
the US poverty hne). Getting back to gini's astrophys- 
ical context, one can consider this question in a more 
visual way by examining the flux distribution function of 
Figure 13. If the distribution of pixels was such that the 
surface brightness threshold cut near the peak of the flux 
distribution function, one would expect the morpholog- 
ical parameters to deviate significantly from the values 
they would have had if the surface brightness threshold 
were lower. If two galaxies at the same redshift have flux 
distributions with identical shape, but different normal- 
ization (i.e., the galaxies have different luminosity), one 
would expect that much less of the fainter galaxy's flux 
distribution would meet the surface brightness criterion 
than the brighter galaxy's. In particular, one would ex- 
pect the the size of the fainter galaxy to be artificially 
small because fewer pixels would be associated with the 
galaxy. Additionally, the gini of the faint galaxy would 
be artificially small because the dynamic range of the 
fluxes in the pixels would be suppressed. In an example 
such as this, the ideal thing to do would be to scale the 
surface brightness threshold so that the threshold would 
apply for the same point on the flux distribution func- 
tion curve, relative to the mean surface brightness, for 
each galaxy. This is the economic equivalent of scaling 
the threshold for calculating gini by the mean income per 
person in each country. 

In order to estimate how much gini and size are affected 
by surface brightness thresholds, we focus on the ends of 
the luminosity distribution. The ten lowest luminosity 
galaxies have a mean surface brightness of {SB)io-u} ~ 
0.015 counts/pix/sec, as calculated using the fiducial 
pixel selection method, mean gini G = 0.332 ± 0.022^, 
and mean size I = 152 ± 18 pixels. The ten high- 
est luminosity galaxies had a mean surface brightness of 

^ The errors listed are the standard deviation of the mean, not 
the sample standard deviation. 
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Fig. 14. — Upper panels: HST / ACS images of {left to right) the low luminosity galaxy BX 1030, the high luminosity galaxy BX 984, 
the low surface brightness galaxy BX 691, and the high surface brightness galaxy MD 169. Each galaxy has typical morphology for the 
galaxies within the sample (i.e., BX 1030 has typical gini and size of the ten lowest luminosity galaxies). The images are stretched on a 
logarithmic scale to match the flux distributions below. Lower panels: Flux distribution functions for the galaxy above. We label the gini 
and size determined with the fiducial pixel selection scheme as "default." In the leftmost two panels, we consider the effects of the surface 
brightness threshold on the gini-luminosity and size-luminosity correlations. The dashed line indicates the surface brightness threshold 
scaled from its default value as the ratio of the mean surface brightness of the 10 most luminous galaxies to the 10 least luminous galaxies. 
Both gini and size are recalculated for the typical high luminosity galaxy BX 984 using the scaled surface brightness threshold, and are 
labeled "L." In the rightmost two panels, we consider the effects of the surface brightness threshold on the gini-surface brightness and 
size-surface brightness correlations. The dashed line indicates the surface brightness threshold scaled from its default value as the ratio of 
the mean surface brightness of the 10 highest surface brightness galaxies to the 10 lowest surface brightness galaxies. Both gini and size 
are recalculated for the typical high surface brightness galaxy MD 169 using the scaled surface brightness threshold, and are labeled "SB." 



{SB)high — 0.03 counts/pix/sec, mean gini 0.527±0.025 
and mean size / = 447±92 pixels, as calculated using the 
fiducial pixel selection method. Representative galaxies 
(low luminosity: BX 1030, high luminosity: BX 984) for 
each sample are shown in Figure 14. Since the mean 
surface brightness of the highest luminosity galaxies is a 
factor of two greater than the mean surface brightness, 
we scale the surface brightness threshold for the high lu- 
minosity galaxies by that same factor of two in order 
to estimate the effects of truncation bias on the gini- 
luminosity and size-luminosity relations. When we select 
pixels for the high luminosity sample using this higher 
surface brightness threshold, we find that the mean gini 
reduces to G = 0.467 ± 0.024, and the mean size to 
I = 292 ± 65 pixels (see the lower panels of Figure 14 
for an example of how the surface brightness threshold 
affects size and gini). Both the mean gini and size of 
the high luminosity sample, as calculated with the scaled 
surface brightness threshold, are > la from the means of 
the low luminosity sample (calculated using the default 
pixel selection method). This suggests that the correla- 
tion between size or gini and luminosity is real, and not 
the result of truncation bias. However, we caution that 
the mean surface brightnesses were calculated using only 
pixels that were above the surface brightness threshold, 
and hence the calculated mean surface brightnesses are 
subject to truncation bias. 

We note that there is also a very tight correlation 
between the gini coefficient and mean surface bright- 
ness (Figure 15), while there is only a weak (2cr) anti- 
correlation between size and surface brightness. In or- 



der to determine if low surface brightness galaxies have 
artificially low ginis due to truncation bias, we raise the 
surface brightness threshold of the highest surface bright- 
ness galaxies by the ratio (a factor of four) of the mean 
surface brightness of the ten highest surface brightness 
galaxies to the mean surface brightness of the ten lowest 
surface brightness galaxies (Figure 14). We then recal- 
culate the ginis of the highest surface brightness galax- 
ies using pixels which meet the new surface brightness 
threshold, just as we did with luminosity. The mean gini 
of the high surface brightness sample (using the fiducial 
pixel selection method) is 0.599 ± 0.024, and the mean 
gini of the low surface brightness sample is 0.227 ±0.009. 
We find that mean gini of the highest surface brightness 
galaxies, calculated with the increased surface brightness 
threshold, decreases by ~ 0.2 to 0.444 ± 0.026, which is 
not enough to destroy any of the gini-surface brightness 
correlation. 

The above trends are consistent with the results from 
simulations of our morphological analysis procedure. As 
described in Section 3.2, in these simulations we added a 
series of artificial object images to the real ACS image. 
In this case, the artificial images consisted of the pixel 
light distributions for objects among the brightest in our 
sample, with the highest gini vales. We then recovered 
the corresponding morphological parameters for simu- 
lated objects in relatively clean locations. Next, we ran 
the same simulations after scaling down the simulated 
object images by factors of 2 and 4 (roughly equivalent 
to lifting the surface brightness threshold by the same 
amount). We found that gini and size for the scaled. 
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Fig. 15. — Surface brightness as a function of gini (left); surface brightness as a function of UV absolute magnitude (right). The small 
black dots denote members of the photometric BX/MD sample. Spectroscopically confirmed BX/MD galaxies are marked with larger blue 
triangles, with members of the spike outlined with blue circles. Red asterisks mark the DRG sample. 



fainter object images are lower by amounts consistent 
with what is outlined above. 

Closer examination of the flux distributions demon- 
strate (c.f., Figure 14) that the flux distribution functions 
of low gini and high gini objects have fundamentally dif- 
ferent shapes. In particular, the most luminous galaxies 
and those galaxies with the highest surface brightness 
have most of the flux concentrated in just a few pixels. 
This means that both the gini coefhcient and average 
surface brightness must necessarily be quite high for high 
luminosity galaxies. 

We also consider the fact that we do not K-correct the 
UV luminosities to a single waveband. Therefore, the UV 
luminosity we use in our analysis effectively samples the 
galactic SEDs at a redshift-dependent wavelength. This 
could affect a correlation in the following sense. If there 
is a tilt to the typical SED, the G band (from which we 
derive UV luminosities) will probe a different part of the 
SED as a function of redshift. Therefore, lower redshift 
galaxies could appear systematically brighter or fainter 
than their higher redshift counterparts. If the morpholo- 
gies had a redshift dependence, this would imprint it- 
self in the luminosity-morphology relations. In order to 
see if underlying size-luminosity and gini-luminosity cor- 
relations exist if the UV luminosity samples the same 
part of the galaxy SED, we perform Spearman corre- 
lation tests on the spike galaxies alone. The central 
wavelengths for the spike range from 1442 to 1455A, a 
range much narrower than that of the entire spectro- 
scopic BX/MD sample. Since the morphologies and lu- 
minosities of spike galaxies are not significantly different 
from those of field galaxies, we expect any morphology- 
related correlation present in the spike galaxies to also be 
present in the larger population. We find that, for the 
fiducial pixel selection scheme, size and UV luminosity 
are still correlated at the 2.7a level, and gini and UV lu- 
minosity are correlated at the 2.5a level. Therefore, the 



UV luminosity-size and UV luminosity-gini correlations 
appear to be real. 

5.4. Comparison with Recent Work 

It is useful to compare our findings to those of Law 
et al. (2007), whose sample of spectroscopically con- 
firmed z - 2 BX/MD galaxies in the GOODS-N has a 
redshift distribution that significantly overlaps with our 
sample. We find that, when we calculate the morpholog- 
ical parameters of our spectroscopic BX/MD sample in 
the same way as Law et al. (2007, "Scheme B"), the mean 
values of size, gini, and multiplicity are almost identical 
to those of the z ~ 2 sample in the GOODS-N field. 

In general, our results are consistent with those of Law 
et al. (2007). In particular. Law et al. (2007) also find a 
strong correlation between gini and luminosity, and size 
and luminosity. Our results on the strength of the cor- 
relation between morphology and stellar age and mass 
are mildly inconsistent with Law et al. (2007), who find 
no correlation between stellar mass and age and mor- 
phology. We find weak evidence for correlations between 
stellar age and mass and morphology; however, the cor- 
relations only marginally significant at the 1.5 — 2fT level. 
The main difference between the two analyses lies in the 
relationship between extinction and morphology. Law 
et al. (2007) find a 3.3a correlation between extinction 
and size, and a 2.4(t correlation between extinction and 
gini. We do not find a correlation between extinction 
and any morphological parameter. However, the method 
of estimating extinction with SED fitting is subject to 
large errors. Therefore, the lack of consistency in the 
strength of the extinction- morphology relationship in the 
two analyses is not meaningful. 

The most interesting result from this part of our anal- 
ysis is that there is no strong correlation between the 
physical properties or star formation history and galaxy 
morphology. This is in stark contrast to what is ob- 
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TABLE 5 

Mean and 1-cr Standard Errors of Morphological Parameters for Each Galaxy Sample 





All 




(SB) < 0.01 count/pix/sec 


0.01 < (SB) < 0.02 


count/pix/sec 


Parameter 


BX/MD 


DRG 


BX/MD 


DRG 


BX/MD 


DRG 


Number of Galaxies 


282 


43 


39 


19 


185 


23 


Size[pix] 


251 ± 10 


215 ± 31 


272 ± 23 


158 ± 31 


252 ± 12 


269 ± 49 


Gini 


0.402 ±0.006 


0.308 ±0.012 


0.269 ±0.006 


0.240 ± 0.010 


0.381 ±0.005 ( 


D.356± 0.010 


Multiplicity 


6.10 ±0.29 


9.36 ± 1.05 


10.44 ±0.86 


12.10 ± 1.81 


6.17 ±0.32 


7.49 ± 1.06 


Flux[counts/sec] 


4.12 ±0.23 


2.51 ±0.41 


2.46 ±0.21 


1.41 ±0.29 


3.54 ±0.16 


3.49 ± 0.67 


(SB) [10-2 counts/pix/sec] 


1.64 ±0.11 


1.17 ±0.36 


0.90 ±0.11 


0.89 ± 0.25 


1.40 ± 0.09 


1.30 ± 0.34 



served at the present epoch, where star formation history 
and morphology are closely tied (Kauffmann et al. 2004) , 
though consistent with the results of Law et al. (2007). 

6. MORPHOLOGIES OF OPTICAL- AND IR-SELECTED 

GALAXIES 

In addition to the sample of optically selected BX/MD 
galaxies, we have identified a sample of infrared-selected 
DRGs, most of which are likely to have a similar range 
of redshifts as the BX/MD galaxies (Franx et al. 2003; 
van Dokkum ct al. 2004; Krick et al. 2006). In this sec- 
tion, we will explore morphology as a function of galaxy 
selection method. 

We use two samples for our morphological compari- 
son. The first sample is the set of BX/MDs that lies on 
our four ACS pointings. This sample includes all photo- 
metric BX/MD candidates galaxies (including those with 
spectroscopic confirmation), and excludes obvious fore- 
ground objects brighter than 22 mag in the F814W filter. 
There is a total of 282 objects in this sample. We use 
this expanded sample for two reasons: it is a factor of 
three larger than the spectroscopic sample, and we want 
to control for the fact that the spectroscopic sample has a 
higher mean flux than the full sample of BX/MD galax- 
ies, and may not be representative of the BX/MD sample 
as a whole. The second sample consists of 43 DRGs with 
Kg < 21 mag. Initial studies of DRGs selected by the cri- 
terion of J — Kg > 2.3 (Vega) have found a mean redshift 
of w 2.5 (e.g., van Dokkum et al. 2003; Forster Schreiber 
et al. 2004). However, recent work indicates that a signif- 
icant fraction of bright DRGs may be at z < 2 (Consclicc 
et al. 2006; Papovich et al. 2006). This population ap- 
pears to be a mix of obscured star-forming galaxies and 
old, passively evolving galaxies (Papovich et al. 2006). 

In order to compare the two samples, we determined 
the mean values of the morphological parameters for each 
sample (Table 5), and performed a 1-dimensional KS 
test for each morphological parameter. We find that, 
although the mean object sizes of the two samples are 
similar, the hypothesis that the size distributions of the 
BX/MD and DRG samples are drawn from the same 
parent population is excluded at the 99.5% level. This 
is consistent with the findings of Law et al. (2007). As 
demonstrated in Figure 16, there are far more DRGs with 
sizes of < 100 pixels than BX/MD galaxies. Both the 
mean values of gini, multiplicity, and total fiux and the 
KS tests performed on these parameters indicate that 
the morphology distribution of the DRGs is significantly 
different from that of the BX/MD sample. In particu- 
lar, while 55 of the 282 BX/MD galaxies have G > 0.5, 
none of the DRGs have ginis that high. There are more 
high multiplicity DRGs than BX/MDs. AU of the high 
multiplicity DRGs have low surface brightness. The total 



fluxes of the DRGs tend to be smaller than the BX/MDs. 

We would like to understand why the morphologies of 
the DRGs are so different from those of the BX/MDs. 
One concern is that the DRG sample is significantly 
fainter, both in terms of flux and surface brightness, than 
the main BX/MD sample. Of our 43 DRGs, only 17 are 
brighter than TZ = 25.5, which is the threshold for the 
BX/MD sample. Of those 17 objects, 9 are also clas- 
sified as BX/MDs. In Section 5, we showed that UV 
luminosity is correlated with both size and gini for the 
spectroscopic BX/MD sample. If hmiinosity is also corre- 
lated with morphology for the photometric BX/MDs and 
DRGs, then differences in the morphology distributions 
could be attributed to the fact that the DRG sample is 
fainter than the BX/MD sample. Similarly, in both sam- 
ples, mean surface brightness is correlated with all mor- 
phological parameters (see Figure 15 for the correlation 
with gini), with the possible exception of galaxy size, and 
so it is possible that the differences in the morphology 
distributions can be associated with the different surface 
brightness distributions of the DRGs and BX/MDs. 

However, we are interested in determining if there are 
differences in the two samples that might be due to some- 
thing other than luminosity or surface brightness. In or- 
der to control for the differences in morphology between 
the two samples that might be attributed to differences 
in surface brightness, we create two subsamples with dif- 
ferent mean surface brightness. The criterion in binning 
the galaxies by surface brightness is that the mean sur- 
face brightness of the BX/MD and DRG subsamples be 
approximately the same for each bin. Therefore, we cre- 
ate surface brightness bins: (SB) < 0.01 counts/pix/sec, 
and 0.01 < (SB) < 0.02 counts/pix/sec. We do not cre- 
ate any higher surface brightness bins because there is 
only one DRG with (SB) > 0.02 counts/pix/sec. We 
present the results of our statistical tests in Tables 5 & 
6. 

For the higher surface brightness bin, the KS tests 

demonstrate that the hypothesis the DRG and BX/MD 
morphology distributions for each parameter are drawn 
from the same distribution cannot be ruled out to signif- 
icant confidence. However, for the lower surface bright- 
ness bin, there are some significant differences. Even 
though the mean surface brightness of both samples in 
this bin are virtually identical, the size and flux distri- 
butions are different at the 99.9% level. The DRGs are 
much smaller than the BX/MD galaxies. The KS test on 
the gini distributions is marginal; the null hypthesis can 
only be ruled out at the 90% level. The multiplicity dis- 
tributions do not appear to be different to a significant 
level. 

We would like to understand why the low surface 
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Fig. 16. — Morphology distributions for the BX/MD sample (blue solid lines) and the DRG sample (red hatched lines). 



brightness DRGs are so much smaller and fainter than 
the low surface brightness BX/MDs. One possible ex- 
planation lies in our choice of observed waveband. At 
z ~ 2.3, the F814W filter probes the galaxies in rest- 
frame UV, and hence, is sensitive to recent, relatively 
unobscured star formation activity. If our main inter- 
est is to use morphology to determine the distribution of 
star-forming regions in a galaxy, the rest-frame UV may 
not trace the star formation of a population that consists 
largely of heavily obscured galaxies, such as the DRGs. 
It is possible that the low surface brightness DRGs have 
large star formation regions, but are heavily obscured by 
dust. DRGs with higher surface brightness could simply 
be less obscured. Alternatively, the low surface bright- 
ness population of DRGs might represent those galax- 
ies which have little active star formation. In either 
scenario, these DRGs might appear smaller and fainter 
in the rest-frame UV than their low surface brightness 
BX/MD counterparts. With both Spitzer observations 
and spectroscopic redshifts, it will be possible to per- 
form robust population synthesis modeling for the DRGs. 



Such modeling will enable estimates of star formation 
and extinction, and indicate the nature of these low sur- 
face brightness galaxies. 

There is a possible systematic effect that might bias es- 
timates of all of the morphological parameters. Our orig- 
inal justifications for using a redshift-independent sur- 
face brightness threshold for pixel selection were that the 
redshift distributions of our BX/MD spectroscopic and 
photometric samples ought be to be similar, and that we 
were forced to use a redshift-independent threshold since 
we do not have spectroscopic redshifts for most of the 
DRGs. However, if the redshift distribution of the DRGs 
were different from that of the BX/MDs, as perhaps in- 
dicated by recent work (Conselice et al. 2006; Papovich 
et al. 2006; Reddy et al. 2006), it would be more prudent 
to use a redshift-depcndent surface brightness threshold. 

The main result from this section is that the UV mor- 
phologies of DRGs appear to be substantially different 
from those of the BX/MD star-forming sample. There 
arc two main contributors to the overall differences in 
morphology. The first contributor is related to sample 
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TABLE 6 

KS TEST COMPARISON OF BX/MD AND DRG MORPHOLOGIES. 



Surface Brightness Criterion Size Gini Multiplicity Flux 

all 0.004549^ 0.000021^ 0.013885* 0.000003* 

(SB) < 0.01 counts/pix/sec 0.003550 0.101108 0.193174 0.001701 

0.01 < (SB) < 0.02 counts/pix/sec 0.514717 0.271975 0.251298 0.129829 

^Thc KS probability that the morphological parameters of BX/MD and DRG samples 

arc drawn from the same distribution. 



selection. Our two samples have different mean surface 
brightnesses, and morphology appears to be correlated 
with surface brightness for both samples. However, even 
if one assembled DRG and BX/MD samples of identi- 
cal mean surface brightness distributions, there would 
be morphological differences between the samples, no- 
tably at the low end of the surface brightness distribu- 
tion. At low surface brightness, the DRGs are substan- 
tially smaller and (necessarily) fainter than the BX/MDs. 
Without more data on the DRGs, it is not possible to 
determine the underlying cause of these morphological 
differences. 

7. DISCUSSION 

We have presented the first robust examination of 
galaxy morphological parameters as a function of envi- 
ronment at z > 1.5. In particular, optical and near-IR 
spectroscopic redshifts for our high-redshift galaxy candi- 
dates were a crucial component for determining whether 
or not the galaxies in our sample were associated with 
the significant overdensity at z = 2.300±0.015. The abil- 
ity to locate these continuum-selected galaxies in redsliift 
space sets this work apart from attempts to study galaxy 
properties as a function of environment in high-redshift 
protoclusters surrounding radio galaxies (e.g., Overzier 
et al. 2006; Kajisawa et al. 2006). In these studies, red- 
shift estimates for non-Lya-emitting galaxies were purely 
photometric, with associated uncertainties of 6{z) > 0.5 
(Overzier et al. 2006). Such uncertainties are too large 
for distinguishing between genuine protocluster member- 
ship, and chance projection. 

Previous modehng of the stellar populations of 1.5 < 
z < 2.9 UV-selected galaxies in the Q1700 field indi- 
cated a significantly higher mean stellar mass and age 
for protocluster galaxies at z = 2.3, relative to galaxies 
in less dense environments (Shapley et al. 2005; Steidel 
et al. 2005). The difference in mean age for protocluster 
galaxies is consistent with theoretical expectations that 
galaxy-scale overdensities should collapse earlier when lo- 
cated within large-scale overdensities. Along with the 
difference in average stellar mass and age as a function 
of environment, we find no corresponding difference in 
the distribution of rest-frame UV morphological proper- 
ties, as described by the non-parametric gini and multi- 
plicity coefficients, and the angular size. We now place 
these results in the context of studies of galaxy properties 
as a function of environment spanning from the present 
epoch back to when the Universe was roughly a third 
of its current age. We also highlight promising future 
observational and theoretical directions. 

7.1. The Lack of a Morphology-Density Relation in the 
Q1700 Field 



While robust environmental trends have been demon- 
strated in the nearby universe for both star-formation 
history and structural parameters, we find only a connec- 
tion between environment and star-formation history in 
and around the Q1700 protocluster. There are many fac- 
tors that can account for the lack of a strong connection 
between morphology and environment for the galaxies in 
our sample at z > 1.5. First, it is important to stress 
that, unlike galaxies in the local universe, for which the 
structural parameters and indicators of star-formation 
history (e.g., color, 4000A break, specific star-formation 
rate) are tightly linked (Strateva et al. 2001), UV- 
selected star- forming galaxies at z > 1.5 do not show the 
same connection between structure - at least in the rest- 
frame UV - and star-formation history. The lack of con- 
nection between morphology and star-formation history 
has been demonstrated for spectroscopically-confirmed 
high-redshift galaxies with modeled SEDs in the Q1700 
(this work, Section 5) and GOODS-N (Law et al. 2007) 
fields. If galaxy morphology and star-formation history 
are uncorrelated at z = 2.3, and star-formation history 
is the property most fundamentally connected with over- 
density, then a strong environmental dependence of mor- 
phology should not be detected. The next significant 
point is that all of the morphological analysis in the 
Q1700 field has been performed on HST/ACS F814W 
images, which probe the rest-frame UV for galaxies at 
z > 1.5. A fairer comparison with morphological stud- 
ies at lower redshift, for which morphological parame- 
ters have been estimated from rest-frame optical images, 
will require imaging our full sample in the near-IR at 
high spatial resolution. Third, it has now been widely 
demonstrated that, unlike galaxies at z — 1, most of 
which are characterized by well-defined centers and cir- 
cularly symmetric morphologies that divide clearly into 
disk and bulge- dominated systems, galaxies at z > 1.5 
generically exhibit irregular morphologies whose physi- 
cal interpretation is unclear (Law et al. 2007; Lotz et al. 
2006; Elmegreen et al. 2004a). Since the range of ob- 
served morphologies at z; > 1.5 is largely disjoint with 
that at lower-redshift, and cannot be divided into clear 
subsets of disk and spheroidal systems, at the very least 
caution should be used when searching for morphologi- 
cal differences as a function of environment at high red- 
shift, in a manner analogous to low redshift environmen- 
tal analyses. 

It is also important to place the photometric proper- 
ties of our sample of high-redshift galaxies on the same 
scale as used for galaxies at lower redshift. In the analy- 
sis of the stellar populations of galaxies in the SDSS as a 
function of environment, Kauffmann et al. (2004) exam- 
ine several properties such as optical color, 4000 A break 
strength, and specific star-formation rate. Analysis of 
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galaxies in the z ^ I DEEP2 survey are presented in 
terms of the average rest-franic U — B color as a function 
of galaxy overdensity (Cooper et al. 2006b), and the av- 
erage fraction of "blue" and "red" galaxies as a function 
of overdensity (Cooper et al. 2006a) or of group vs. field 
environment (Gerke et al. 2006). Here, "blue" and "red" 
refer to the two portions of the bimodal distribution of 
galaxy colors, observed both locally (Strateva et al. 2001) 
and at 2 ^ 1 (Cooper et al. 2006b). 

In particular, using the set of best-fit population syn- 
thesis models for UV-selected galaxies in the Q1700 field 
(Shapley et al. 2005), we compute the distributions of 
rest-frame U — B color and 4000 A break strength, 
£>„(4000)^ for spike and non-spike galaxies, respectively. 
The mean rest-frame optical color is [U ~B) = 0.67[0.59] 
for spike [non-spike] galaxies, while the correspond- 
ing mean value of the 4000 A break is (£>„(4000)} = 
1.13[1.20]. Comparison with the color-magnitude distri- 
butions presented in Cooper et al. (2006a) indicates that 
both spike and non-spike galaxies fall squarely within the 
color distribution of "blue" galaxies in the DEEP2 sam- 
ple. The i)„(4000) values for both spike and non-spike 
galaxies are at the low extreme for galaxies in the SDSS, 
significantly below the transition range of ~ 1.6 — 1.7, 
which divides the bimodal distribution of galaxy stellar 
populations at z ~ (Kauffmann et al. 2003). Even 
the oldest galaxies in the Q1700 BX/MD sample, in- 
cluding those contained in the protocluster, would be 
classified as "blue" galaxies using DEEP2 and SDSS cri- 
teria. Furthermore, the stellar population differences be- 
tween spike and non-spike galaxies are not significant 
when compared with the differences between "blue" and 
"red," or early-type and late-type galaxies, in lower red- 
shift surveys. These points are consequences of the young 
age of the universe at z = 2.3, the fact that the sample 
under consideration was selected in the rest-frame UV, 
requiring the presence of current star formation, and the 
fact that, at z ~ 2, unlike in the local universe, galaxies 
with stellar masses of ~ 1O^^M0 still harbor active star- 
formation (Erb et al. 2006; Papovich et al. 2006; Reddy 
et al. 2006). 

We must also take into account the fact that UV- 
selection leads to an incomplete census of the range of 
galaxy stellar populations at z ~ 2 (Reddy et al. 2005; 
van Dokkum et al. 2006), missing galaxies with red 
rest-frame UV colors or faint rest-frame UV luminosi- 
ties stemming from dust extinction or an evolved stellar 
population. It has been shown that the DRG criteria 
yields a sample of high-redshift galaxies largely comple- 
mentary to the set of BX/MD galaxies. While the DRG 
sample appears to be composed of a mixture of dusty 
star- forming, and passive galaxies (Papovich et al. 2006), 
the typical stellar M/L ratio is significantly higher in 
this population than in the BX/MD sample. Ideally, we 
would like to construct complete spectroscopic samples 
of both BX/MD and DRG objects in the Q1700 field, 
and characterize the typical stellar populations and mor- 
phologies of combined BX/MD-I-DRG samples in both 
spike and non-spike environments. The associated limi- 

^ This estimate of the magnitude of the 4000 A break is used 
by, e.g., Kauffmann et al. (2003, 2004), for chara<;terizing stellar 
populations in SDSS, and is defined as the ratio of flux density at 
4000 - 4100A to that at 3850 - 3950A. 



tation is the difficulty of obtaining spectroscopic redshifts 
for DRGs, which are significantly fainter in the rest-frame 
UV than BX/MD objects of the same mass, and char- 
acterized by lower-equivalent-width Ha emission in the 
rest-frame optical, except in the case of an AGN (van 
Dokkum et al. 2004). Current and future multi-object 
near-IR spectrographs on 8 — 10 meter-class telescopes 
will aid in the effort to obtain spectroscopic redshifts for 
DRGs. Without spectroscopic redshifts, it will not be 
possible to characterize the environments of DRGs, due 
to the magnitude of error in typical photometric redshifts 
(van Dokkum et al. 2006). By including passive galax- 
ies in the mix, we will be able to study the fraction of 
galaxies without ongoing star formation as a function of 
environment, which is more analogous to lower-redshift 
studies of the relative fraction of "blue" and "red" galax- 
ies than our previous investigation of "blue" and "slightly 
less bhie" objects. Indeed, we have only shown so far 
that BX/MD galaxies in the spike appear to have started 
forming stars earlier than non-spike galaxies. However, 
since galaxies must have ongoing star formation to sat- 
isfy the BX/MD color criteria, we have little informa- 
tion about the past or future epoch when star formation 
ceases, and the related causes. 

We have presented several explanations for the lack 
of observed dependence of BX/MD morphology on envi- 
ronment, including the fundamental differences between 
low- and high-redshift galaxy morphologies, the appar- 
ent lack of connection between irregular morphology 
and star-formation history for UV-selected star-forming 
galaxies at z ~ 2, and the young stellar populations of all 
the galaxies considered here when compared with sam- 
ples at lower redshift. However, other recent observations 
of morphologies of galaxies at similar look-back times 
suggest evolution in morphology as a function of redshift, 
which may also be relevant to the expected environmen- 
tal dependence of morphology in high-redshift protoclus- 
ters. Indeed, Law et al. (2007) find that the mean value 
of the gini coefficient is significantly higher for a sample of 
z ~ 3 LBGs, relative to that for the corresponding sam- 
ple of UV-selected galaxies at z ~ 2, indicating more nu- 
cleated morphologies for the higher-redshift sample. Es- 
timating typical sizes of star-forming galaxies from z ^ 1 
to z ^ 5, Ferguson et al. (2004) conclude that galaxies at 
higher redshift are smaller, at fixed luminosity. As a re- 
flection of the acceleration of structure formation within 
the protocluster, Steidel et al. (2005) suggest that spike 
galaxies attain a level of maturity at z ^ 2.3, that will 
be reaciied by typic;al "fiekr' (i.e., non-spike) galaxies at 
at slightly lower redshift. Given the physical overdensity 
of the Q1700 protocluster, the morphological properties 
of spike galaxies should be similar to those of non-spike 
galaxies at (z) ^ 1.8, as opposed to non-spike galaxies 
at (z) = 2.3 (as in the current non-spike set). However, 
our small sample size, and the relatively small offset in 
effective redshift (Az ~ 0.5), most likely prevent signifi- 
cant detection of differences in the average gini or size of 
spike galaxies due to their lower effective redshift. 

7.2. The Existence of a Color- Density Relation in the 
Q1700 Field 

There are multiple lines of argument to explain the lack 
of difference in the distributions of rest-frame UV mor- 
phologies for spike and non-spike galaxies in the Q1700 
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field. While there is no correlation between environ- 
ment and morphology for our sample of galaxies, the 
stellar populations of spike galaxies do appear older and 
more massive than those of galaxies in less dense envi- 
ronments. As discussed above, these differences are more 
subtle than the low-redshift distinction between "red" 
and "blue" galaxies, yet the fact remains that we do find 
statistically significant environmental differences in the 
stellar mass and age distributions. It is important to 
consider this result along with those from the DEEP2 
survey, that the average color is redder and fraction of 
red galaxies is higher in group, or denser, environments 
at z < 1.3, but that there is no detectable dependence 
of galaxy color on local environment at z > 1-3 (Cooper 
et al. 2006a; Gerke et al. 2006). Given the ~ 1 Gyr 
timescale for galaxy colors to transform from "blue" to 
"red" following the cessation of star formation, the ap- 
pearance of environmental trends at z = 1.3 indicates 
that the preferential formation of "red" galaxies in denser 
DEEP 2 environments must commence at ^; < 2. For the 
z = 2.3 progenitors of the galaxies probed in the DEEP2 
survey, then, no environmental trends are expected. 

Our results at z > 2 can only be reconciled with the 
DEEP2 results if the galaxies in the z = 2.3 protoclus- 
ter do not represent the progenitors of the bulk of the 
DEEP2 sample, but rather the progenitors of galaxies 
residing in richer environments by z = 1.3. This appears 
to be the case, if indeed the z — 2.3 protocluster viri- 
alizes to become a ~ IO^^Mq cluster by z ~ (Steidel 
et al. 2005). We investigate this question with the aid of 
the Millennium Run simulation (Springel et al. 2005), a 
large, high-resolution N-body simulation of cosmic dark 
matter structure formation. According to dark matter 
halo merger trees constructed from the output of the 
Millennium simulation (Springel et al. 2005; Lemson & 
Virgo Consortium 2006), the typical mass ratio between 
a ^ IO^^Mq dark matter halo at z ^ and its most 
massive progenitor at z = 1.3 is roughly a factor of ~ 6. 
Therefore, by z = 1.3, the Q1700 spike galaxies should 
be residing in a halo of mass ~ 2 x 10"Mq. The DEEP2 
survey is more heavily weighted towards galaxies resid- 
ing in the field or in less massive structures. The bulk 
of groups featured in the DEEP2 sample of Gerke et al. 
(2006) range in mass from 5 x 10^^ < Af < 5 x 10^'^, with 
a negligible fraction of systems with M > 10^^ Mq. The 
DEEP2 sample of Cooper et al. (2006b, a) probes a simi- 
lar range of environments. Therefore, the environment of 
the likely z = 1.3 descendant of the z = 2.3 protocluster 
is not well-represented within the DEEP2 c^iivironmcntal 
samples, and there is no contradiction between our ob- 
served color dependence on environment at z = 2.3 and 
the lack thereof in the DEEP2 survey at z = 1.3. 

7.3. Future Directions 

There are several promising theoretical and observa- 
tional directions to pursue. More work is required to 

determine the detailed three-dimensional distribution of 
the protocluster galaxies, and the manner in which they 
are grouped. Given that the angular distribution of 
galaxies contained within the spike in redshift space ex- 
tends over at least ^ 11 comoving Mpc on the sky, and is 



not obviously segregated from non-spike galaxies (Stei- 
del et al. 2005), the protocluster defined in redshift space 
does not constitute a single, well-defined, virialized ob- 
ject. Therefore, we would like to understand the mass 
of typical structures containing protocluster galaxies at 
z = 2.3 that will inhabit ~ lO^^M© clusters by z ~ 0. 
Based on the Millennium simulation merger trees, it ap- 
pears that the z ~ 2 progenitor halo mass distribution 
of ^ 1O"^^M0 clusters at z ^ typically contain between 
one and five ^ IO^^Mq dark matter halos. Since the dark 
matter halo mass in which a galaxy resides affects the 
amount and fate of mass and gas it accretes (Keres et al. 
2005; Dekel & Birnboim 2006), a careful understanding 
of the way in which protocluster galaxies populate dark 
matter halos at z ^ 2 will shed light on their future 
star-formation histories. While detailed dark matter plus 
hydrodynamics simulations of massive cluster formation 
(Kravtsov et al. 2005) have been performed, most of the 
close comparisons with observations have been tuned to 
z ~ 0. We would also like to employ a close comparison 
of the detailed star-formation histories and structures of 
simulated protocluster and field galaxies prior to z = 2.3. 

Furthermore, future morphological studies of statisti- 
cal samples of z ~ 2 in the rest-frame optical will provide 
a cleaner view of galaxy structure. The camera currently 
on HST with the minimum acceptable angular resolution 
for this type of study is NIC2 on the Near Infrared Cam- 
era and Multi-Object Spectrometer (NICMOS), though 
its small field of view provides a survey efficiency too low 
for building up a statistical sample of rest-frame opti- 
cal morphologies. However, the Wide Field Camera 3 
(WFC3), to be installed in HST during Servicing Mis- 
sion 4, will enable this type of study in future years. Our 
detailed observations of the environment of the Q1700 
protocluster, plus a close comparison with simulations of 
cluster formation tuned to z = 2.3, will provide a unique 
opportunity to learn about the seeds of the environmen- 
tal dependence of galaxy properties. 
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